On hitting times of affine boundaries by reflecting 
Brownian motion and Bessel processes 



Firstly, we compute the distribution function for the hitting time 
of a linear time-dependent boundary t t— > a + bt, a > 0, b G R, by 
a reflecting Brownian motion. The main tool hereby is Doob's for- 
mula which gives the probability that Brownian motion started inside 
a wedge does not hit this wedge. Other key ingredients are the time 
inversion property of Brownian motion and the time reversal property 
of diffusion bridges. Secondly, this methodology can also be applied 
for the three dimensional Bessel process. Thirdly, we consider Bessel 
bridges from to with dimension parameter S > and show that 
the probability that such a Bessel bridge crosses an affine boundary is 
equal to the probability that this Bessel bridge stays below some fixed 
value. 
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1 Introduction 

Let B = {Bt : t > 0} denote a real valued Brownian motion and / : 
[0, +oo) i— >• R, a "nice" function. Introduce the first hitting time of / by B 
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via 

H f := inf{i : Bt = f(t)}. 

There are many important applications (e.g. in sequential analysis, see, e.g., 
Lerche [21], and mathematical finance, especially in the theory of American 
options, see, e.g., Peskir and Shiryaev [Mj Chapter VII) which call for the 
distribution of Ht for some particular functions /. For linear boundaries 
fit) = a + bt we refer to Doob |13j and Durbin [T3], for square root bound- 
aries f(t) = \J a + bt see Breiman [8], Ciffarelli et al. [10], Shepp [33], Delong 
[H],[l2] and Yor [37], and for parabolic boundaries fit) = a + bt 2 , Groene- 
boom [19] , Salminen [29] and Martin-L6f [22] . We refer to Durbin [15] (with 
appendix by Williams) for a practical method via integral equations to calcu- 
late the hitting probabilities of time dependent boundaries. For the method 
of images, see Lerche [21j . In Alili and Patie [3] a new approach based on 
functional transformations of processes is developed and its connection with 
the method of images is studied. 

In this paper we compute the distribution of the hitting time of a linear 
boundary for 

• reflecting Brownian motion (RBM), 

• three dimensional Bessel process (BES). 

Our approach is based on a formula by Doob presented in |13] , as well as 
the time inversion property for Brownian motion and Bessel processes, and 
the time reversal property of diffusion bridges. In fact, the case of RBM has 
been already studied via Doob's formula and the time inversion by Abundo 
([2]; Corollary 3.4) for straight lines with positive slope. For Bessel processes 
with dimension parameter 5 > results for hitting times of straight lines 
through origin can be found in Pitman and Yor ([25]; Section 8 p. 332) 
and in Alili and Patie ([!]; Theorem 5.1) - we review these results below in 
Remarks 1151 and Theorem 1161 respectively. 

The above description shows that one of the aims of our paper is to 
relate formulae about these (classical) first hitting times which appear to be 
rather scattered in the literature. Moreover, seemingly different formulae 
for the same densities may be found in different papers but this discrepancy 
can often be explained, e.g., via the Poisson summation formula. 

We now present our main notation. All processes we study are defined 
on the canonical space C of continuous functions uj : R+ i— > R. Let 

Ct := (j{oj s := w(s) : s < t} 
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denote the smallest u-algebra making the coordinate mappings up to time t 
measurable and take C to be the smallest cr-algebra including all cr-algebras 
Ct, t > 0. Processes are identified in (C,C) via their probability measures; 
e.g., P BM , P BBM , and P BES ^ denote the probabilities under which the 
coordinate process with u = x is a Brownian motion, a reflecting Brownian 
motion and a Bessel process with dimension parameter 5 , respectively. 
Throughout the paper the notation 

H a}b := inf{t : oj t = a + bt} 

is used for the hitting time of a linear time dependent boundary. Notice that 
for x = 0, by the scaling property of Brownian motion and Bessel processes, 



d \b- 2 H abA , b>0, 

^" 2 fla|6|,-l, b<0, 

and if f a ^ denotes the density of H a ^ then 



'b- 2 f ab ,i(b 2 u), b>0, 
b- 2 f a \b\,-i{b 2 u), b<0. 



fa,b(u) = \^_ 2 , ,, 2 , L ^ n (2) 



Also for x = 0, the time inversion property - a basic tool in our approach - 
is valid for Brownian motion and Bessel processes with dimension parameter 
5 > initiated at (for < 5 < 2 the Bessel process is assumed to be 
reflecting at 0), that is, under Pq M and p BES ^ 

{co t ■ t > 0} = {tu 1/t :t>0}. (3) 

We remark that this property can also be used to derive the distribution of 
the last hitting time from the distribution of the first hitting time (and vice 
versa). Indeed, for a, b > let 

G bta := sup{t : oj t = b + at}. 
Then, under P^ M and P^ ES{5) , it holds 



sup{t : ojt = b + at} = sup{t > : toj-y/ t = b + at} 

b 

= supjt : u 1/t = - + a} 
= 1/ infj-u : oj u = a + bu}, 
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and, hence, 

G b , a = 1/H a>b , (4) 
in other words, e.g., for a, b > 

P° ES(5) (G b , a < t) = Pf^> [H a , h > i) . (5) 

Next we recall well known formulas for the hitting times of constant 
boundaries by RBM and BES. Firstly, 

E* BM (exp(-A# O)0 )) = ( cosh(aV2A)) ~* (6) 



from which using series expansion and inverting term by term we obtain 

oo 

P RBM (H a ,o e du) /da = -£= J2 (4* " 1) e"^- 1 ) 2 ^. (7) 

fc=— oo 

Making an elementary parity manipulation one can find a different but 
equivalent form for the density of H a< o 

oo 

P? BM (#« o G du) /du = -=L= V (-l) fc (2A; + 1) e -(2fc+D 2 a 2 /2« (8) 

V27rn 3 fc= _ oo 

displayed in [7] p. 355 and 641. 

For three dimensional Bessel process it holds 

(exp(-A^o)) = -r^L, 0) 
smh(av 2AJ 



and inverting 

P* ES ®(H afi edu)/du (10) 



K = — OO 



^ ((2fc + l)V - it) e~ 



(2fc+l) 2 a 2 /2u 



which coincides with the formula presented in [7] p. 463 letting therein x \, 0. 
We refer also to Pitman and Yor [27] where the subordinators, denoted by 
Ct and St, generated by the positive powers of the right hand sides of © 
and ©, respectively, are studied. 

For ease of the reader (and of ourselves !), when comparing certain theta 
function type series indexed by Z or N, we indulge in stating the following 
elementary lemma. 
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Lemma 1. Let tp be a function defined on odd integers 2k + 1, k G Z. 
Assume that there exists e such that tp(2k + 1) = e ip{—2k — 1) for all k G Z. 
Then 

oo oo 

^ ^(2n + 1) = (1 + e) J>(2fc + 1), 

n=— oo k=0 

where it is assumed that the sum on the RHS is finite. In particular, if ip is 
even, i.e., e = +1 then 

oo oo 

^(2n + l) = 2j>(2fc + l), 
n=— oo fc=0 

and if ip is odd, i.e., e = — 1 i/ten 

oo 

^(2n + l) = 0. 

n=— oo 

The outline of the paper is as follows: Doob's formula is discussed in 
the next section, and in the third section the distribution of H a f> f° r RBM 
is presented. In the fourth section the BES case is treated, and the paper 
concludes with a fifth section, where the distribution of the maximum of 
a general Bessel bridge is connected with the distribution of H a i, for this 
Bessel bridge. 



2 Doob's formula 

In his famous paper on the Kolmogorov-Smirnov test, Doob [13J derived the 
expression for the probability that Brownian motion started inside a space- 
time wedge does not hit this wedge. We first give some symmetry properties 
of the function characterizing this probability. 

Proposition 2. For a, (3, a, b > let 

G{a,(3;a,b) := P^ M {-(a + /3t) <uj t <a + bt V t > 0) . (11) 
Then G has the symmetry properties 

G(a, p; a, b) = G(/3, a; 6, a) = G(a, b; a, /3), (12) 
and for all c> the scaling property 

ct a 

G(a,P;a,b) =G(-,f3c;~,bc). (13) 
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Proof. The first equality in (|12p follows from time inversion: 



G(a, /3; a, b) 



jyBM 

r 



(a + /3t) <uj t <a + bt V t > 0) 
(a + fit) < tujy t <a + bt V t > 0) 



pr (-(f +/3)<«i/i 

G(/3, a; 6, a). 



< - + b V f > 



The second equality in (|12p is obtained by the spatial symmetry of Brownian 
motion. For (1131) consider 



G(a, P; a, b) = (-(a + pt) < uj t < a + bt V t > 0) 



^0 
pBM 



(a + /3c 2 t) < uj c 2 t < a + 6c 2 i V t > 0) 

-(- + Pet) < - oj c 2 t <- + bct V t > 
c c c 



= G(-,/?c;~ 6c), 
c c 

where Brownian scaling has been used. 

We now state Doob's formula; for the proof, see |13j . 
Theorem 3. For a, P, a, b > 

oo 

G(a, P;a,b) = l-J2 (e" 2A + e~ 2Bfc - e" 2Cfc - 



□ 



-2D 



k ) , (14) 



k=l 



where 



A k := A fc (a, P; a, b) := /e 2 &a + (k — I) 2 Pa + fc(A; - l){ba + a/3), 



B k (a, P; a, b) := (k - l) 2 ba + k 2 Pa + k(k - l)(ba + a/3), 
C k (a, P; a, b) := k 2 (ba + /3a) + k(k - l)ba + k(k + l)a/3, 
L> fc (a, /3; a, 6) := & 2 (>a + Pa) + fc(fc + l)6a + k(k - l)a/3. 



Remark 4. T/ie symmetry properties of G can, of course, be seen also from 
the expressions for A k , B k ,C k , and D k . Indeed, 

A k (a, P; a, b) = B k (a, b; a, 0), 

and, since, 

C k (a,P;a,b) = k 2 (ba + /3a) + k 2 (ba + a/3) + k(aP - ba) 
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and 

D k (a, j3; a, b) = k 2 (ba + flat) + k 2 (ba + a/3) — k{aj3 — bat), 
it is easily seen also that 

C k (a,(3; a, b) = D k (a, b; a, /3). 

Moreover, 

A k (a, (3; a, b) = A k ((3, a; b, a), B k (a, f3; a, b) = B k (f3, a; b, a) 

and 

C k (a,(3; a, b) = C k (f3, a; b, a),D k (a, (3; a, b) = D k ((3, a; b, a). 

We also note that each of A k ,B k , C k and D k satisfies the same scaling prop- 
erty as G. Notice also that if we introduce 

T(m, n) := m 2 ba + n 2 j3a + mn(ba + a/3) 

then 

A k = T(k,k-l) and B k = T(k-l,k). 

The limiting distribution function (found by Kolmogorov) of the Kolmogorov- 
Smirnov statistics is 

oo 

F K {\) : =l + 2^(-l) fc e- 2fc2A2 , A > 0, 

k=l 

and was identified by Doob to be the distribution function of the maximum of 
the absolute value of a standard Brownian bridge. Moreover, Doob observed 
that this distribution function is closely related to the probability that RBM 
never hits a straight line with positive slope. Indeed, taking in (llip a = a 
and /3 = b gives 

A k = B k = (2k - l) 2 ba and C k = D k = (2k) 2 ba, 

and yields this probability as displayed in the following 

Corollary 5. For a,b > it holds 

P^ 7 (\u t \<a + bt V t > 0) = P* BM (H atb = oo) 

oo 

= l + 2^(-l) fc e- 2fc2ab (15) 

k=l 

= e*{2ab/Tr), 



7 



where 

6» :=yV-l)*e- 



kez 

is the "modified" Jacobi theta function. 
Remark 6. Trivially, the LHS of [T5\) equals 

P^ M (supflwtl -bt)<a 
\t>o 

so that a \— > @*(2a 2 /it) appears as a distribution function and equals Fx- 
For a more general discussion for such suprema involving Bessel processes, 
see Theorem\HK We note that the scaling property of Brownian motion also 
yields that the LHS of t!5\) is a function of the product ab only. 

Using time inversion Abundo [2] (see also Scheike [32]) computed the 
probability that a Brownian bridge stays inside a space-time wedge. Let 
x, y € R and u > be given, and denote by the probability measure 

associated with Brownian bridge from x to y of length u. The next result is 
given in [2], Proposition 3.5 and Theorem 3.3. To make the exposition more 
self contained we give here a proof which differs slightly from the proof in 

Theorem 7. For u > and a, (3, a, b > 

,BM 



F o,u,y H a + Pt)<u>t<a + bt V t e (0, u)) 



Gla^ + ^a^+^l, (16) 



where — (a + (3u) < y < a + bu and G is as defined in ( |j^[ ). In particular, 
for a = a, (5 = b, and \y\ < a + bu 

Po B ^(M<« + ^ Vt€(0,u)) 

oo 

= V (-l) n exp(-2n 2 a(6+-) + 2na^ . (17) 
^— ' V ii ill 



n=— oo 



Proof. Using 

(Wf, o<t< U }, p™ y ) I (to + hi o < t < u}, pgft) 



s 



we write 

,BM 



Vo,% v (-a-Pt<u>t<a + bt Vte(0,u)) 



p o,u,o - fa < at + -* < a + bt V t € (0, «) 
p £*fo - C9 + |)t < wt < o + (6 - V i € (0, «) 



To proceed recall 

(K;0 < i < u},Pg5) = N ^w«t/(«-t);0 < * < ^} 5 P? M ) 
and, hence, 

p o!^o - C9 + < <*t < a + (P - ~)t V t e (0, u 

= P BA/ f -a - 09 + < — u ut/M < a + (6 - V t E (0, u)) 

= P**(-a-(0+ V^<— <*<a + (&- V,>0 

\ M V + U V + U U V + U 

= P? M (-a - (a + fiu + y)- <u v < a + (a + bu-y)- V u > 
\ u u 

n( a , a + y h , a ~y 
= G [a,p H ; a, 6 H 



by the symmetry property of G. This proves formula (|16p and (|17p follows 
immediately. □ 

From (JTTJ) , integrating over y leads to the distribution function oiH a b, a,b > 
0, for RBM, see Abundo [2] Corollary 3.4. In the next section, see Theorem 
[9l we find this distribution also for negative values of b. 



Remark 8. The fact that variables a, b, y and u appear on the RHS of |P. 
in such a particular way can be explained again (cf. the last sentence in 
Remark^ via the scaling property. Indeed, from the proof above (consider 
the case u = 1 ) 

Po!K,(NI<a + ta Vie (0,1)) 

= P^ M {-a - (a + b + y)v < u v < a + (a + b - y)v V«>0) 
= P^ M (-a - (a + b + y)v < cco v/c 2 <a + {a + b-y)v V v > 0) 

= P? M (-- - (a + b + y)cs <u s <- + {a + b- y)cs V s > o) . 
\ c c / 
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Choosing c = a yields 

>BM 

0,1,2/ 



Pl¥A\u t \<a + bt Vt 6(0,1)) 



P^ M (\oj s + ays| < 1 + a(a + 6)s V s > 0) 



showing that the LHS of |7?j) can be viewed as a function of b' := a(a + b) 
and \x := ay only. 

3 Distribution of H a ^ for a reflecting Brownian 
motion 

In the next theorem we give an explicit form of the distribution function of 
H a)b , a > 0, b G R, for RBM initiated at 0. For b > formula dTS} can be 
found in Abundo [2] Corollary 3.4 (however, there seems to be a misprint in 
formula (3.6) in [2]; in the exponential n should be replaced by n 2 ). Clearly, 
if RBM is initiated from above the line then the distribution of H a ^ can be 
obtained from the distribution of H a ^ for BM. 

Theorem 9. For a > and b € R and u > suc/i that u < —(a/b) for 
b < and i/iere is no constraint if b > it holds 

V* BM (H a>b > u) 

£ ( _ 1} V 2 ^ ( $ ( Q + 2ka \-^f- a + b ^ 2ka )), (18) 



oo 



A;=— oo 

where $ denotes the standard normal distribution function, i.e., 

$(*):= f -Le^dy. (19) 

J-oo V 2-7T 

JTie density, with the same constraints, is given by 
P* BM (H a , b G du) /dn 

-. +oo 

=■ A a , fe (n). (20) 
Moreover, in case b > 

p?-'K..>^(r" ( ^ >2?) - < 2i > 

1, ify<a, 
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and for y > a 

Pq™ E dy) /dy = l^aA^y, (22) 

in case b < 

Pf-K.„<v)4f» KBM( ^ > ^ ) ' (23) 

1, itu> a, 



and for y < a 



^ BM ["B a , b E dy) /dy = |^|A a , 6 (^). (24) 



Proof. We need to prove the result for b < (for 6 > see [2]). For this, by 
scaling, it is enough to consider iJ 0) -lj i-e., we take b = — 1 (cf. (P)). Using 
formula (|16|) in Theorem [7] we have for u < a 



P BBM (Ha,.! >u)= T U P BBM (ff 0j _! > u , Wu G dy) 
J 

= f " Pq™ (#„,_! > n | w u = y) P BBM (oo u G dy) 

JO 

= f " P??„ ("(« " s) < u s < a - s Vs G (0,u)) 

j — (a— «) 

x P BM (ui u G dy) 



P^o (-(« - - «)) < "S < s - (u - a) Vs G (0, u)) 

(a— u) 

x P bm (oj u G dy), (25) 

where, in the last step, we have used the time reversal property of diffusion 
bridges, see Salminen [30J , i.e., for Brownian bridge we have 

({UH :0<t< u}, V™ y ) = ({uj u - t :0<t< u},P B ^ x ) . (26) 

Using spatial homogeneity and (|16p in Theorem [7] we obtain 

Ku,o (-(« -(u-a))<u 8 <s-(u-a) Vs G (0, u)) 

= P??-» -{u-a))-y<u s <s-{u-a)-y Vs G (0,u)) . 

= p o!«f-y (-(« - u + y + s ) < u s < a - u - y + s Vs G (0, u)) . 

= G(a -u + y,-;a-u-y,-). 

u u 
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To write down the explicit form of the function G we use (|14[) , and to stress 
the dependence of A k , B k , and D k on y we use notations A k (y), B k (y), 
Cfc(y), and D k (y), respectively. It holds 

A k (y) = - (k 2 (a - u - y) + (k - l) 2 (a - u + y) + 2k(k - l)(o - u)) 
u 

= -{2k - lfa + - ((2k - l) 2 a 2 - (2k - l)ay) , 

B k (y) = ^[(k~ lf(a - u - y) + k 2 (a -u + y) + 2k(k - l)(o - u)) 
= -(2k- ifa + - ((2k - 1) V + (2k - l)ay) , 

C k (v) = - (2k 2 (a -u) + k(k-l)(a-u-y) + k(k + l)(o - u + y)) 
u 

= -(2k) 2 a + - ((2k) 2 a 2 + 2kay) , 



and 

D k (y) = - (2k 2 (a - u) + k(k + l)(a - u - y) + k(k -l)(a-u + y)) 
u 



i2k) 2 a + - ((2k) 2 a 2 - 2kay) 



Now we have, see (j25j) . 

P$ BM (H a ^ > u) (27) 
[ 1 - E (e" 2Afe(s/) + e- 2Bftfe) - e" 2 ^ - e~ 2D ^) J ^(y) dy, 



-(a—u) 

where 



fc=i 

1 



Mv) :=-^=^ l 2u . 



Straightforward integration yields: 
f " e- 2A ^ Vu (y)dy 

J ~{a—u) 

f a ~ u 1 / v 2 2 \ 

= / -== exp -f ((2fc - l) 2 a 2 - (2k - l)ay) dy 

J-(a-u) V2vru V 2u u J 



,2i2/.-l)-„ / 1 e -(»-2(2fe-l)a) 2 /2« rfy 



-(a- u ) \/2vru 

e 2(2fc-i)2a ^ ^ a-^-2(2fc-l)a ^ _/-(«-«)- 2(2A: - l)a 



v 7 ^ 



=: A(fc,it), 
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e- 2B ^My)dy 

— (a—u) 

2(2fc-l) 2 a ( a ~ u + 2 ( 2k ~ * f~( a ~ u ) + 2 ( 2k ~ l)a 



(2fc-i) 2 a f$ ^ a-u + 2(2k- l)o \ _ $ 
: B(fe,«; 



(a—u) 

(a — u) + 4A:a 



=: C(fc,u), 



and 



(o— u) 

_ p 2(2fe) 2 a ^ ( a-u-Aka \ f -(a-u) - 4ka 

= : D(fc,u). 
Hence, we have 

P* BM (#a,-i > u) 

oo 

= 1(0, u) - (A(k, u) + B(k, u) - C(k, u) - T>(k, u)) , (28) 

k=l 

where 

T/n \ f a ~ U \ ^( Ol-U 

1(0, u) = $ — =- - $ 



Rearranging the terms in (|28|) yields (|18p with 6 = — 1, and differentiating 
term by term leads to the density given in (|20p . □ 

Remark 10. Taking in U8\) 6 = and differentiating yields formula §3$ 
for the density of the hitting time of a point for RBM. Clearly, the density 
given in $20\) also converges as b — > to (EJ). 
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4 Distribution of H a ^ for a 3-dimensional Bessel 
process 

In this section we derive an expression for the distribution function of H a ^ 
for a 3-dimensional Bessel process initiated at 0. The probability measures 
associated with a Bessel bridge from x to y of length u with dimension 
parameter 3 and a killed (at 0) Brownian bridge are denoted by «,f an d 
P^^ J , respectively. Recall that, in fact, 

TiKBM _ r,BES (9 q\ 
*■ x,u,y r x,u,y 

which follows from BES being Doob's /i-transform of KBM with h{x) = x, 
i.e., for T t G C% and y > 

(see McKean [23], Williams [36], Revuz and Yor [28] p. 450 and Borodin 
and Salminen [7] p. 75). However, for some arguments below it is good to 
have a different notation since these bridges are conditioned from different 
processes. 

We start with a proposition which says that a Brownian bridge from 
x > to y > conditioned not to hit is identical in law with the killed 
Brownian bridge. 

Proposition 11. For x,y,u > 

P™ = P ^(-K>° V0<t<u). (31) 
Proof. For 0<ui<---<u„<u consider 
P SU„ Ki G dzx, ...,u Un €dz n ;u t >0 V < t < u) 

= Pf M G cfei, . . . , w Un € dz n , w n £dy; H >u) /Pf M (w« G dy) 
= Pf BM (w Ul G dzi, . . . , w Un G dz n , w u G dy) /Pf M (oj u G dy) , 

where Hq := inf{t : oJt = 0}. Similarly, 

P2?„ ("t > V0<Kk) = Pf BM K e dy) /Pf M K G dy) . 
Consequently, 

G dzi, . . . ,u Un G d^ n \uh > V < t < u) 
= P« BM (oj Ui G dz u . . . ,oj Un G d^,^ G dy) /Pf BM K G dy) 
= P Sf K G dz\ , . . • , w n „ G dz n ) . 
Then, the monotone class theorem implies the stated identity. □ 
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Theorem 12. For a > and b G R suc/i i/iai u < — (a/b) for b < (an<i 
i/iere is no constraint if b > 0) it holds 

= 2 [ e - (• (^^) - * (-7=)) 



— a — bu ( a + bu — 2ka 



u \ \/U 



(32) 



where (p is the standard normal density. 



Proof. Consider first the case b < and u < —(a/b). Without loss of gen- 
erality we may take b = — 1. Proceeding as in the proof of Theorem [9] (cf. 
(j25j) ) we have 



ra— u 

P$ ES >u)= P** S Q (io s < s + a - u Vs G (0, u)) 

JO 

x P£ BS K G dy). (33) 

By (J29J) and @TJ, 

P^K< S + a-«Vse(0, U )) 



PoT K < « + a " « G (0,u)) 



P^(0<w,<» + a-uV«e(0,«)) 

= lim — ^— : „, . ; ; — . 34 

X V> P w K > Vs G (0, u)) V ' 

From Theorem [7] formula f)16|) we have for x,y > 

Pjjs (0 < uj s < s + a - u Vs G (0, «)) 

= PoSf*- w H<w s <s + a-«-j/VsG(0, u)) 

I X CL — X \ 

= G[y,-,a-u-y, , (35) 

\ u u J 



and 



P^K>0VsG(0,u)) 
\BM 



lim P^ 1 (-y < ^ < JV - y Vs G (0, u)) 
lim G<( y ,-;iV- y ,^y (36) 
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From (|36p we may recover the well known formula (cf. [7] p. 174) 

Vy,u,x K > Vs G (0, u)) = 1 - e~ 2xy / u . (37) 

To write down the explicit form of the function G in (j35j) we use (j!4[) . The 
functions A^, B^, and are now denoted by A^x, y), B/ t (x, y),Ck(x, y) 
and Dk(x,y) and it holds 

^4fc(x, y) = - [k 2 (a - u - y)(a - x) + (k - ifxy 
u 

+k(k - l)((a — u — y)x + (a - x)y)] , 



Bk{x, y) = - \{k -l) 2 (a-u- y)(a - x) + k 2 xy 

+k(k — 1) ((a — u — y)x + (a — cc)y) ] , 

Ck(x, y) = - [k 2 ((a-u- y){a - x) + xy) 
u 

+k(k — l)(a — u — y)x + k{k + l)(a — x)y\ , 

and 

D k {x, y) = - [k 2 ((a-u- y)(a - x) + xy) 
u 

+k(k + l)(a — u — y)x + k(k — l)(a — . 



Moreover, 



A fe (0, y) = D k (0, y) = - (fc 2 (a - it) - , 
u 

B fc+ i(0, y) = C fe (0, y) = - (k 2 (a - u) + %) 



and 



^:_^ fe ,)_(*-l)£Zi + 2, 
ax u u 

tk- g x (x,y)--D k .-- ^ (x,y)--k ^ . 
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We calculate the limit in (fM|) using l'Hospital's rule: 



P y,u S o (us < s + a - u Vfl€(0,lt)) 

1- Er=i (e" 2Afc + e~ 2Bfc - e~ 2C * - e~ 2D *) 

~ xio i - e - 2x y/ u 

_ B[e-^ + g=i (K - D'k) e- 2 ^ (0) + (g±i ~ C' k ) e~^) 

y/u 

2(a - u) 

^ A I 



^ ^ Z^a ^ ^ e 2k 2 a e -2k 2 a 2 /u ^ Q -2kay/u _ Q 2kay/u 

=1 

<x> 

+£ 



^2k 2 a^—2k 2 a 2 /u t^—2kay/u _|_ ^2kay/u 



k=l 

oo 



2(a — u) ^^2k 2 a e -2k 2 a 2 /u e -2kay/u 
oo 

+ £ 



^ k=—oo 



e 2k 2 a e -2k 2 a 2 /u e ~2kay/u 



k=—oo 

u 2kay / u 



A _ 2fc(a - m) \ e 2fc2 ae _ 2fc 2 a 2 / 



fc=— oo 



where, in the second step, -Bi(O) := Bi(0,y) = 0,^(0) := Afc(0, y), and 
Cfc(0) := Ck(0,y). Now we can proceed by evaluating the integral given in 
(133]) . Recall that 

rr S ("u£dy) = -^=e-y 2 ^dy, 
and consider the integrals 

h := / e -2k 2 a 2 /u e 2kay/u _J^ = & -y 2 t 2u dy 
JO V 27TU 3 



^ p -fa-2fca) 2 /2 M| 

o v^vra 3 



and 



y \/27rn 3 



a—u 



2k{a-u) I -^ =Q -{y^af/2u dy 

lO V 27TM 3 
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Writing 



/— 2y(y-2ka) ^ 2ka?/2u ^ 

7o 



+2/ca / 

JO 



2 y e -(y~2ka) 2 /2u 



V2-KU 3 



dy 



--■ hi +/12, 



and 



I 2 = 2/ca 



2 ?/ r -fa-2fca) 2 /2^ 



\/2~KV? 

■ hi — I22 



dy 



-2k [ 
Jo 



^ e -(y-2fca) 2 /2« 

\/2iru 



dy 



it is seen that I\2 = hi- Straightforward calculations yield 



a — u ( a — u — 2ka 



'11 



u 



■LP 



+ 2 $ 



a — u — 2ka 



2ka 
fa 



Furthermore, 

JO V27TU 

+2(2fc) 2 a 



2.. I " ^ =e -(y2kaf/2u dy 
V27TU 

a — u — 2/ca 



+2(2&f a $ 



a — u — 2ka 



2ka 
fa 



Consequently, 

00 

P^ ES (H a ,-i > u) = £ e 2fc2a (h ~ h) 



k=—oo 
00 



/ 2fcu — a + u 



fc=— 00 

00 



„2fc 2 a 



2 £ 

fc=— c 

+ (1 + (2£;) 2 a) ^ 
18 



95 



a — u — 2ka 



a — u — 2ka 



2/ca 
7^ 



which is the claimed formula for b = — 1, and, by scaling (see ([I])) we may 
extend this for general 6 < 0. 
Let now 6 = 1, and consider 



P$ES ( HaA >u )= f pBES ^ s<s + a y s£ (o, u )) 
Jo 

x P^ ES (co u e dy). (39) 

We have (see (fMj) ) 

P?S? (w.<a + oVa€(0,«)) 



,. Pji(0<w,<* + aV»6(0 )U )) 

hm D , f ; ; • (40) 

xio P™(w s > Vs e (0,u)) 



Using (JT6J) gives (see (|B3|)) 



Pfi (0 < w s < s + a Vs € (0,«)) 

= P o,u,y-x (~ x < u s < s + a - x Vs <E (0, u)) 

G(*,^;a-x,l + ^ 
it u 

To find the explicit form of the function G, let Ak(x,y), Bf c (x,y),Ck(x,y), 
and Dk(x,y) denote now the terms A^, Bk,Ck, and in (fT4"|h It holds 

y) = - [fc 2 (a + u - y)(a - x) + (k - l) 2 xy 



u 



+k(k - l)((a + u — y)x + (a - x)y)] . 



Bk(x, y) = - \{k -l) 2 (a + u- y)(a - x) + k 2 xy 
u L 



+k(k — 1) ((a + u — y)x + (a — x)y) ] , 



y) = - [k 2 ((a + u- y)(a - x) + xy) 
it 



+/c(fc - l)(o - x)y + k(k + l)(a + u- y)x] 



and 



D k (x, y) = - \k 2 ((a + u- y)(a - x) + xy) 
u L 



+k(k + l)(a — x)y + /c(/c — l)(a + u — y)x] 
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Notice that, e.g., A k (x, y) can be obtained from A k (x, y) by substituting — u 
instead of u and multiplying by —1. 
Moreover, 



A k (0,y) = C k {Q,y) = £ (k\a + u) - ky) , 
4+1 (0, y) = £> k (0, y) = - (k 2 (a + u) + ky) , 



and 



^( X y)-- k ^±^ + l 

ax u u 

dB k a + u y 
—— = (k - 1 h -, 

ox u u 



dC k , dD k .a + u 

x ^y) = — a — ( x ' y) = k — • 



dx ' dx 
Consequently, as in (f38|h 



C^ <fl + i Vt€(0,«)) 



-2k 2 a ^—2k 2 a 2 / u f^—2kay/u _ ^2kay/u 



k=—oo 

and, further, 



1+ 2(Q + U) ^> 

OO 

_|_ e -2k 2 a e -2k 2 a 2 /u ^2kay/u _|_ g -2A <//; ,; 1 
fc=l 

^ / _ 2fc(q + n) ^ e _ 2 fc 2 a e -2fc 2 a 2 /« e 2fea ? ;/ tt; ^-Q 



fc=l 

oo 



P (^ ail >n)-2^e ^ -j= j 



a + u — 2ka\ _ / 2/ca\ 



+ 1 - (2*0^) $ = - $ 



- 



and, by scaling, see ([T]), we may extend this for general 6 > 0. Since formulas 
(|2ip -(|24 p follow easily, the proof of the theorem is now complete. □ 



Remark 13. It is easily checked that taking b = in llSty) and differentiating 
yields formula $10\). Moreover, letting y \, in fi41\ ) and (38\) . respectively, 
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we obtain the following formula with b = ±1 (if b < then the formula is 
valid for u < a/ (—b) ) 

P^ S (oj t < a + bt Vt6(0,u)) 

= 1 + 2 E f 1 " 4fc 2 fl(a + H ) e ~^a(^u)/u^ (42) 

k=l ^ u J 

The case with general b can be proved, e.g., via scaling. Notice that this 
probability is a function ofa(a + bu)/u only. In fact, see Theorem [F, 



+ « Vie (0, „)) = P*** sup a* < J^±M 

\t<l » U J 

Putting in fi42\ ) 6 = gives the well known formula for the distribution of the 
maximum of a Brownian excursion found in Chung f^, see also Kennedy 
\20y . Durrett and Iglehart 176]/ . Biane and Yor fdjj and Pitman and Yor \26^ . 
The distribution of the maximum of a general Bessel process is discussed in 
Theorem 17 

Next we consider the distribution of H a<b for a three dimensional Bessel 
process initiated above the line 1 1— > a + bt. In this case, the density can be 
found using (local) absolute continuity of the probability measure induced 
by the Bessel process with respect to the probability measure induced by a 
killed Brownian motion. 

Theorem 14. For x > a > and b > it holds 

)BES 



(a + bt)(x - a) c _ (a _ x _ bt) 2 /2t ^ 

In particular, 



(H a b £dt)= K ' ™l e -{a-x-btr/2t dt (43) 

txV2M 



?E ES (Ho, b € dt) = ~^f hHl2 dt. (44) 

Proof. Since H ab < Hq for b > 0, we obtain, using absolute continuity ([30]) . 

P* ES (H a , b E dt) = ; H a , b e dt 

= £±*P«« (H a , b E dt) 
x 



{a + bt)(x-a) ^ (a _ x _ bt) 2 /2t 
tx\j2irt 
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by the well known formula for the density of the hitting time of straight lines 
for Brownian motion (see, e.g., [7] p. 295). For formula (j44p . recall that 
is an entrance boundary point for the three dimensional Bessel process and, 
hence, (|44"|) is obtained by taking in (|4"5]l a = and letting x — > 0. □ 



Remark 15. For Bessel processes with dimension parameter 5 > 2 we have 
that Hq i, > P BES -a.s. for all b > (see Shiga and Watanabe \3J$ Theorem 

3.3). On the other hand, it is also well known that uit/t — > P x -a.s. 
for all x > 0. Consequently, < H a ^ < oo a.s. for all a > and b > 0. 
From Pitman and Yor ]25$ Section 8 p. 332 the following result for hitting 
times of lines going through the origin holds 

P s x {H , b edt)/dt = (bt) u e-^ b2t+x2/t )/{2tx' / K u (xb)), x>0, (45) 
and 

P 5 (#o,& G dt) /dt = (j) T-V^/rH, (46) 

where v = (5 — 2)/2 and K v is the modified Bessel function of the second 
kind (see Abramowitz and Stegun p. 374). Since Ki/ 2 (x) = \fne~ x / 'y2x 
and r(l/2) = t/tt it is seen that formulas |^5| j and fl^ffi ) coincide with |^3[ ) 
and ^44\ ) i n case ^ = 3 and a = 0. 

Alili and Patie derive in [4j an expression for the density of hitting time 
of a Bessel process to a straight line. We wish to compare their formulas 
in case of RBM and BES with the formulas presented in Theorems [9] and 
[T2l The next result is Theorem 5.1 in [3] (we have, however, corrected a 
misprint in [4] by putting in minus signs in front of b and b 2 in the exponent 
of the first exponential). 

Theorem 16. For a Bessel process of dimension parameter 5 > and initial 
state < x < a it holds 

i( a 3_ a J ) _^ t 

*> BES ^{H a ^dt)/dt = f 



(i+ b -ty+ 2 

where b £ R, v = {5 — 2)/2 and is the Bessel function of the first kind 
and 2v,k is the ordered sequence of the positive zeros of J u . In particular, for 
RBM starting from 0, i.e., v = — 1/2 and 

J-\/2(z) = \ — cos z, Ju 2 (z) = \ — sinz, 

' \ TTZ 1 \ TTZ 
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it holds 



P$ BM (H a , b £dt)/dt 



s/Eia + btfl 2 

oo , 

x^(-l) fe+1 (fc--)7r exp (-((*-£)* 

k=l ^ 



2 7 7 2a(a + 6i), 
= : A,,&(t). (48) 

Comparing the results in Theorems [16] and [U we obtain, with our nota- 
tions in (jH]) and ([20]) . that 

£> a , 6 (i) = Aa, 6 (t). (49) 

This is a Poisson summation type identity to which we would like to devote 
attention in a forthcoming work [31j. For the moment, we only refer to 
Feller [18] and Biane, Pitman and Yor [5] for some related discussion. 



Figure 1 about here. 



Figure 2 about here. 



5 Results for general Bessel processes 

In this section we identify first the probability that a Bessel bridge crosses an 
affine boundary with the probability that the maximum of the Bessel bridge 
is below some value. For brevity, we change the notation used above and let 
Pa denote the probability of a Bessel process with dimension parameter 
5 initiated from x > 0. Moreover, Pq 5 ^ denotes the probability measure 
of a Bessel bridge from to of length u. For < 5 < 2 the Bessel 
process is assumed to be instantly reflecting at 0. Secondly, we show that 
the probability for a Bessel process to stay below the straight line 1 1— )■ a + 1 
equals to the probability that the maximum of the corresponding Bessel 
bridge from to of length 1 is less than y/a. 
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Figure 1: The density of the distribution of the first hitting time oft \— > 3 + i 
by RBM started from 0. 

Theorem 17. For a > 0, u > 0, and b < a/u 

F o,l,o ( su pH + W < a )= P c£,o ( supH} < - • (50) 

\t<« / \t<u / 

Proof. Consider first the case 6 > and let A; denote the left hand side of 
(IKUjl . We have 

A, = P^ (co t <a-bt V0<t<u) 

>(5) 



Po% )0 ( w * < ° - hu + 6i V0 < t < 

,(-5 




Pn 5) {oJ t <a-bu + bt V0 <t <u\u u = 0) , 
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A 3> -i(f) 




Figure 2: The density of the distribution of the first hitting time oft^2> — t 
by RBM started from 0. 



where in the second step the time reversal property of diffusion bridges is 
used (cf. (gSD ). 

By the time inversion of Bessel processes we obtain 

A, = P[, 5) (tuj 1/t <a-bu + bt V <t <u\ uu 1/u = 0) , 

= P { S) L 1/t <{a-bu)- t +b V ~ > ~| u 1/v = OJ , 

= PjP (ou s < (a-bu)(s + ~) + b V s >0 
= (u s <(a- bu)s + - V s > 
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Letting A r denote the right hand side of (|50p we write 

A r = P^i j0 (u t < y/a{a - bu) V0<t<u) 



t &i/t — ■>/ a ( a ~ bu) V — > — | uui\/ u = J 

< ^/ a(a — 6u) s V s > — | Wjy^ = 0^ 

< - 6u) ( s + ^) Vs>oj 

< Va(a - bu) (s + ^\ V s > j . 



= Pjp \^uj s < ^a(a - bu) \ s + V s > Oj . ( 5 1 ) 

Recall that Bessel processes enjoy the Brownian scaling property: 
for all c > under 

{sfcu t/c : t > 0} = {oo t : t > 0}. (52) 
Applying this with c = (a — bu)/a in ([ST]) yields 

A r = P { Q S) (u> as /(a-bu) <as + ^ \f s > 
= Po 5) (ut < (a - + - V t > 
which equals A/. 

For 6 < the proof is slightly simpler since we do not use time reversal. 
Indeed, letting again A; and A r denote the right and left hand side of ([50]) . 
respectively, we have by time inversion 

A i = P c?l (wt < a - bt V < t < u) 



= 



P? } ( w s < as + VOO). (53) 



11, 



For A r it holds (again) by time inversion 

A r = p{P < v / a(a - bu) ( s + ^) V s > j . 
Using hereby scaling (]52p with c = a /{a — bu) we obtain 

A r = p[, 5) (w( a _ 6u ) s / < (a - bu)s + - — — V s > 
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and substituting t = (a — bu)s/a it is seen that A r = A/, thus completing 



Remark 18. We note that for b G R, the result of the theorem may be 
stated as 



where Mf, := sup t<1 {o; t — bt}. 

The next result is a slight extension of Exercise 3.10 in Revuz and Yor 
|28j where reflecting Brownian motion case is considered. We present here 
also the proof for the readers' convenience. 

Theorem 19. For Bessel processes with 5 > it holds 



In other words, letting R = {R t : t > 0} denote a Bessel process started 
from with dimension parameter 5 > and R ' 1,0 = {R% '° : < u < 1} 
the corresponding Bessel bridge from to of length 1. Then 



the proof. 



□ 





(54) 




(55) 



Proof. Consider for a > 




P (Vt, R\2 t < a + X 2 t) 
P (Vt, XR t < a + \ 2 t) 



by scaling for any A > 0. Choosing A 2 



a yields 



P ( sup{R t - t} < a ) = P Nt, y/aRt < a + at) 
\t>o J 




i.e. 




(56) 
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We now note that 



su p|rV7 1 = supU 1 ~ u ) R u/(i-u)}- 
t>o L 1 + t J u <i 



= sup^' 1 ' }, (57) 

U<1 

where, in the first step, we use the substitution u = t/(l + i), and, for the 
second step, we use the well known representation 

{R°/'° : < u < 1} = {(1 - u) i? w /(l- u ) : < u < 1}, 

which is inherited from the corresponding Brownian bridge representation. 
Claim (1551) follows now from (1561) and (1571). □ 



Example 20. For 3- dimensional Bessel bridge the distribution of the supre- 
mum is known. In fact, it was proved by Vervaat I35f that this supremum 
is identical in law with the range of a standard Brownian bridge, i.e., 



sup<^' \ = supB t > -mfB t ' 

t<u L > t<u *<" 



The distribution of the range can be deduced from Feller [11], and we have 
(seefflp. 174) 

.2" 



p( s ^-g"^<v)-i + *E(»-^) 

Consequently, from Theorem\T^ for d = 3 

oo 

P (sup{R t - t} < a) = 1 + 2 V (1 - 4fc 2 a) e" 2fc2a . 
^° ^ 

Example 21. For 1-dimensional Bessel bridge, i.e., reflecting Brownian 
bridge, the distribution of the supremum is also known. From [T^ p. 333 
formula 1.1.8 we obtain 

\ oo 

? 0,«,0| \ = ! _,_ 9 ^ f e -2(2fc)V/« _ e -2(2*-l)V/u" 



pfsup| J B t ^ |<^=l + 2^ 
anrf, hence, from Theorem [TP] /or 5=1 



Po(sup{|Bt| - t} < a) = 1 + 2 V (V 2 ^ - e^" 1 ) 2 ") 

oo 

= l + 2^(-l) 



oo 

\k e ~2k 2 a 



k=l 
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This formula may also be obtained from [T8\) by letting u — > oo and coincides 
with < T75)) . 

Acknowledgement. We thank Mikhail Stepanov for some numerical 
computations and creating the plots of A a fe, see Figure 1 and 2. 
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